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k*" ' Abstract 

. 

\Q , In this paper, we construct a generalized Darboux transformation for nonlinear Schrodinger 

equation. The associated TV-fold Darboux transformation is given both in terms of a summation 

^ \ formula and in terms of determinants. As applications, we obtain compact representations for 

the iV-th order rogue wave solutions of the focusing nonlinear Schrodinger equation and Hirota 

— _ 

equation. In particular, the dynamics of the general third order rogue wave is discussed and shown 
to exhibit interesting structure. 
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1 Introduction 

Darboux transformation, originated from the work of Darboux in 1882 for the Sturm-Liouville equa- 
tion, is a powerful method to construct solutions for integrable systems. The theory is presented by 
several monographs and review papers (see [61 [191 ED- ^ n literature, various approaches have been 
proposed to find a Darboux transformation for a given equation, for instance, the operator factoriza- 
tion method [T], the gauge transformation method [20 |. 114 } IT], the loop group transformation [22] and 
so on. 

It is remarked that Darboux transformation is very efficient to construct soliton solutions. Indeed, 
through iterations, one is often leaded to compact representations in terms of special determinants 
such as Wronskian or Grammian for N— soliton solutions. Such N— soliton solutions are appealing 
both form the theoretical viewpoint and form practical application viewpoint as well. 
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In addition to the soliton solutions, rational solutions are also interesting and Darboux trans- 
formation may be adopted for this purpose. For the celebrated Koteweg-de Vries (KdV) equation, 
Matveev [IB] introduces the so-called generalized Darboux transformation and the positon solutions 
are calculated. Recently, rogue waves, appeared in oceans, have been studied and applied extensively 
in other fields such as Bose-Einstein condensates, optics and superfluids and so on. The very first 
model for rogue waves is the focusing nonlinear Schrodinger (NLS) equation 



which has been an important integrable equation. While the simplest rogue wave solutions are cal- 
culated by Akhemidiev and coworkers, and the construction of higher order analogue is one of the 
challenges as remarked in [5]. In this regard, as pointed out by Dubard et al [8], the solutions, ob- 
tained by Eleonskii, Krichever and Kulagin [10], represent a class of multi-rogue wave solutions. It is 
remarked that the construction method proposed in [8j [9] is very specific and technical, so may not 
be easy to apply to other models. 

The aim of the paper is to propose a simple method to construct multi-rogue wave solutions. 
The main tool is the generalized Darboux transformation. We will reexam the Matveev's generalized 
Darboux transformation for the KdV equation and derive it in a way that could be easily extended 
to other models. Then, we apply the idea to the focusing NLS equation and work out a formula to 
generate multi-rogue wave solutions. 

This paper is organized as following: In section two, we propose a new way to derive the generalized 
Darboux transformation for KdV equation. In the section three, we first apply the proposed method 
to the NLS equation case and obtain the corresponding generalized Darboux transformation for it. 
Then, we reformulate the iV-th fold generalized Darboux transformation in terms of determinants. 
Also, we provide formulae for the iV-th order rogue wave solutions for NLS equation and Hirota 
equation. With the help of this formulae, we consider the dynamics of the general third order rogue 
wave. The spatial-temporal pattern of the solution could form as a triangle or pentagon. Final section 
contains some discussions. 

2 Generalized Darboux transformation for KdV 

Let us first recall the well known classical Darboux transformation for the KdV equation. Consider 
the Sturm-Liouville equation 




(1) 



(2) 



and introduce the following first order operator 



T[i] = d x 
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where *$>i is the fixed solution of (|2|) with A = Ai. Then Darboux transformation 
converts the equation (J2]) into 

-*[lU + u[l]*[l] = A*[l] (3) 

where 

u[l] =u-2QnV 1 ) xx . 

and Wr(^i, = ^i^ x — ^i tX *& is the standard Wronskian determinant. 

The most interesting point here is that one could iterate the above Darboux transformation. 
Indeed, the iV-times iterated or iV-fold Darboux transformation yields Crum theorem 

-^ xx [N] + u[NMN} = X^[N], u [iV]=«-2(lnWr(*i,... 

where 

1 1 Wr(*l,-" ,* N ) ' 
and V&i, ■ ■ ■ , \l/ 7v are solutions of ([2]) at A = Ai, • • • , Ajv, respectively. 

It is obvious that ^i[l] = T[l]^i = 0, namely \Pi is mapped to trivial solution. This fact 
implies that a seed solution may not be used more than once when considering the iterations for 
Darboux transformation. However, as pointed out by Matveev and Salle |19j. a generalized Darboux 
transformation does exist. Let us derive this result in a way which may be readily generalized. To 
precede, we assume that 

*2 = *i(fci + e), 

where k\ = f{\\) is a monotonic function and e is a small parameter. Expanding ^ m a series in e 

^ 2 = ^ 1 + ^ ] e + ^ ] e 2 + ■■■ , 
where W± = ^ 9 g^i^ |fc=fci- Since ^[l] = ^ill] 1 ^ is a special solution for (J3j), so is ^slll. Taking limit 



for this solution, we find 



* ll i]=i™ Tmi(hl+e) =mM>, 



which is a non-trivial solution for (J3j) at A = Ai- This solution may be adopted to do the second step 
Darboux transformation, that is 

T[2] =d x - %4r> «[2] = " 2(lnWr(*i, ^ ] )) xx . 
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Combining these two Darboux transformations, we obtain 

-* M [2] + u[2]tf [2] = A#[2], [2] = g^Cgii^iig) . 

Wr(*i,^^ 1J ) 

This process may be continued and results in the so called generalized Darboux transformation 
for the system ([2]). Indeed, let ^i, • • • , be n different solutions for (|2j) at Ai, A2, ■ • • , A n , and 
consider the expansions 

%(h + e) = m i {k i ) + ^ ] e + --- + ^ ] e m * + --- , k i = f(X l ) (i = 1, 2, • • • , n) 

then we have the following 

Proposition 1 f7&| / 



with 



Wo 

u[N] =u- 2(ln(W 1 )) xx , V[N] = — ± 



solve 



-■* xx [N]+u[N]®[N} =A*[n], 
where m\ + m2 + • • • + m n = TV — n, rrij > 0, rrij £ Z. 

The generalized Darboux transformation presented above may be used to generate both solitions 
and rational solutions for the KdV equation. Let us illustrate this with the following examples. It is 
well known that, the KdV equation 

u t - 6uu x + u xxx = 0, (4) 

takes ([2]) as its spatial part of the spectral problem and the corresponding temporal part reads as 

*t = -m xxx + 6u^ x +3^*. 

In the case of N distinct spectral parameters, we will have the Wronskian representation for the 
N— soliton solution. To get rational solutions, one starts with the the seed solution u = c, c is a real 
constant, and ^1 = sin[A;i + {Akl + 6c)t) + P(h)}, ki = VAT^, P(h) is a polynomial of k\. Now 
expanding the function V&i at k\ = and taking P{k\) = for convenience, we have 



W 1 = (x + 6ct)ki + 



--{x + Qctf + 4t 
6 



Y^(x + Qctf - 2t(x + Qctf 



k\ + 
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therefore 



^f l =x + 6ct, ^ = — (x + 6ctf + 4t, * 



i 



[2] = 



— {x + Qctf -2t(x + 6ct) 2 



then the generalized Darboux transformation provides us the following rational solution for the KdV 



G = 12[279936c 5 t 5 (216c 5 t 5 + 360c 4 xi 4 + 270c 3 x 2 t 3 + 120c 2 x 3 i 2 + 35cx 4 i + 7x 5 ) 
+38880c 4 t 4 (180t 2 + 7x 6 ) + 25920c 3 xi 3 (x 6 + 180t 2 ) + 1620c 2 x 2 t 2 (x 6 + 720£ 2 ) 
+60ci(x 9 + 216(tt 2 x 3 + 4320t 3 ) + 43200i 3 x + x 10 + 5400x 4 t 2 ], 

H = 3888c 4 t 4 (12c 2 t 2 + 12cxt + 5x 2 ) + 4320c 3 t 3 (x 3 + 3t) 

+540c 2 xi 2 (x 3 + I2t) + 36cx 2 t(x 3 + 30t) + x 6 - 720t 2 + 60x 3 £. 



*i = sin[A;i(x + (Ak\ + 6c)t) + P(ifei)] = <if J fci + tf^A;? + *f J A;f + • • • + *f J fe 2Af+1 + • • • . 



In particular, the positon solutions for the KdV equation may be found this way and for a detailed 
analysis we refer to the work of Matveev [17] . 

Concluding this section, we mention that, for those equations which possess solutions represented 
in terms of Wronskians, there are alternative manners to explain the limit solutions [19 | I15 |. [T3] . 

3 Generalized Darboux transformation for NLS 

In this section, we extend the idea discussed above section to the NLS equation and construct a gener- 
alized Darboux transformation for it. Furthermore, we shall show that such Darboux transformation 
enables one to obtain, apart from the soliton solutions, rational solutions including multi-rogue wave 
solutions. 

The focusing NLS equation, equation ([I]) is the compatibilty condition of the following linear 
spectral problems 



equation (jH) 



«[3] = C + T7o 



where 



More general solutions of rational type may be obtained if we expand 



[iC<ri + iQ]tt, 
iCVi + iCQ + ^i(Q x 



iQ 2 ) 



(5b) 



(5a) 



where 
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3.1 Generalized Darboux transformations 

The Darboux transformation in this case is defined as (see [B] and the references there) 



[1] = T[l]tf , q[l] = q + 2(d* - Ci)(P[l]U, 
where ^ 

r[i]=c-cr+(cr-ci)Pi[i], m = ^ w 

and is a special solution of the linear system ([5a|) - ([5bj) at ( = Ci> (-Pl[1])21 represents the entry of 
matrix Pi [1] of the second row and first column and "f " denotes the matrix transpose and complex 
conjugation. 

If N distinct seed solutions ^^(fc = 1,2, ...,N) of (|5ap -(|5b |) are given, the basic Darboux trans- 
formation may be iterated. To do the second step of transformation, we employ Vl/ 2 which is mapped 
to # 2 [1] = T[1]| C=C2 # 2 - Therefore, 

*[2] = r[2]*[l], q[2] = q[l] + 2(C 2 - C 2 )(P[2]) 2 i, 

where 

T[2] = C-C| + (CI-C 2 )P[2], P[2] = ||j^i|. 
In general case, we may have 

Theorem 1 Let . ^9 . • • • . ^at beN distinct solutions of the spectral problem (I5ap - (|5bp at ^ . • • • ,Cn 
respectively, then the N-fold Darboux transformation for NLS equation ^ is 

N 

*[N]=T[N]T[N-1]-..T[1]% q[N]=q[0]+2^2(C:-Ci)(P[i\)2i (7) 



i=l 



with 



9,[i - 1] = (T[i - l]T[i - 2] ■ ■ ■ r[l]) | ( . ( ,*„ ,[0] = 

We remark that the iV-fold Darboux transformation given by ([7J) is equivalent to the determinant 
representation presented in [p9] as will be seen in the next subsection. 

Now we manage to find a generalized Darboux transformation. As in last section, suppose that 
$2 = ^i(Ci + (5) is a special solution for system, then after transformation we have ^[l] — Ti[l]^ r 2- 
Expanding ^2 at Ci> we have 

^(d + 6) = + + vfjV + • • • + ^[ N] 5 N + • • • , (8) 
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where *f ] = ^^*i(C)lc=Cx- 

Through the following limit process 

lim \m]k t+ sm = lim h+t.pii^,]*, = jipn.fi = , lW (9) 

<5^0 5 <5->0 6 

we find a solution to the linear system (|5aj ) -(|5b |) with q[l] and C = Ci- This allows us to go to the next 
step of Darboux transformation, namely 

7\[2] = C - Ci + (Ci " Ci)^i[2], q[2] = q[l] + 2(£ - Cl)(^l[2]) 2 i, (10) 

where 

gi[l]gi[l]t 
nH *i[l]t^[l]' 

Similarly, the limit 

lim [^ + r 1 [2](Ci)][ < y + r 1 [i](Ci)]^ 2 = 

<5->0 5 2 

*! + [TxtlKCi) + T 1 [2](Ci)] +T 1 [2](Ci)T 1 [l](Ci)^ 21 = *i[2], 

provides us a nontrivial solution for the linear spectral problem with q = q[2] and C = Ci- Thus we 
may do the third step iteration of Darboux transformation, which is following 

Ti[3] = c-Ci* + (cr-ci)iM3], a [3] = ^pjf^pr ^ 3 ] = ^ 2 ]+ 2 (cr-ci)(A[3]) 2 i. (n) 

Continuing above process and combining all the Darboux transformation, a generalized Darboux 
transformation is constructed. We summarize our finding as 

Theorem 2 Let \l/i(Ci)) ^2(C 2 ), "j ^n(Cn) be n distinct solutions of the linear spectral problem 
(EaD-dSH) and 



*i(Ci + S) = *i + #f ] <5 + ttJV + • • • + tff 1 * 1 ^ + • • • , (fc = 1,2, • ■ ■ , n) 

6e i/ieir expansions, where 

r i 1 

*• = 7T^J V "' Kj - - (* = 1, 2, ■ ■ ■ ,n, j = 1, 2, ■ ■ ■ )■ 

Defining 

T = r n r n _i---rir , r i = r i [Tni]---r i [i], (*> 1), r = i, (12) 
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where 

[<* + Ti[k] c=Q ] ■■■[5 + Ti[2] c=Ci \[6 + r i [i]| c=fi ]iv. 1 (Ci + 5) ■ ■ • ri(Ci + *)r *i(Ci + s) 



H>i [k] = lim 



<5->0 <5 fc 



*i + E E ( r * ■ • ■ r * [<] • • • ^ ■ ■ • r« [4?]) Ic=c^ 



s=1 mi>h^>->h^>l, 
i>gi>—>gi>l 
if gi = i, then h^ 1 ' < k 

(1 < k < rrii), then the following transformations 

n rrii n 

*[N]=TV, g[iV] = g + 2EE(C-Ci)(^K])2i, (N = n + ^m k ) (13) 

i=l j'=l fc=l 

constitute a generalized Darboux transformation for the NLS equation. 

We remark here that the solution formulae (J7J) and (|13[) . represented in terms of summations, have 
certain merit. Indeed, for non-zero ^k,k = 1, 2, ■•• , N, all the denominators of P[i] and .Pi [7] are 
easily seen to be non-zero in these forms, therefore both (J7J) and (|13p supply nonsingular solutions. 
The former could lead to iV-soliton solutions, while the latter may yield rouge wave solutions. 

Let us consider an example to illustrate the application of above formulae to construction of higher 
rogue wave solutions. To this end, we start with the seed solution q$ = e 1 *. The corresponding solution 
for the linear spectral problem at £ = ih is 

vM^l^-^-y-*) (14) 

l (C 2 e A - C ie - A )ei' J ^ ' 

where 

(h-VhP^T) 1 / 2 „ (h + Vh^l) 1 / 2 A r- , 

Let h = 1 + f 2 , expanding the vector function ^i(f) at / = 0, we have 

^ 1 (/) = vI/ 1 (0) + ^ 11 / 2 + --- (15) 



where 

*i(0) 



-2t + 2ix-i) e-5 1 
(2it + 2x + 1) el 1 * 



[i] - |t + \ - 2ta 2 + fx 3 + ft 3 - 2kc£ 2 - ix 2 + 2tx + it 2 ] e - ^* 

\ [|sc + f t - \ - 2\tx 2 + fx 3 - ft 3 - 2xt 2 + x 2 + 2ixi - t 2 ] e\ 



if 
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It is clear that ^i(O) is a solution for (f5a]) - ([5b|) at ( = i. By means of the formula ([9]), we obtain 



*i(0)*i(0)t' 



^[1] = Hrn ^Ul^^il = ^[1]^ +m(0), 



T ^ = 2i V * l( o)t* l( o);- 

Substituting above data into (fTUl) yields the second-order rogue wave solution 

G l +itG 2 



9[2] 
where 



1 + 



H 



Gi = 36 - 288x 2 - 192x 4 - 1152tV - 864t 2 - 960t 4 , 
G 2 = 360 + 576x 2 - 192t 2 - 384x 4 - 768x 2 t 2 - 384t 4 , 

H = 64t 6 + 192i 4 z 2 + 432t 4 + 396t 2 + 192i 2 x 4 - 288i 2 x 2 + 9 + 108x 2 + 64x 6 + 48x 4 . 

which is nothing but the solution first constructed by Akhmediev et al [5]. The higher order rouge 
wave solutions may be calculated similarly thus we have a general approach to produce these solutions. 

3.2 Determinant forms and higher-order rogue wave 

In generic cases, iterated Darboux transformations may be given compactly by means of determinants 
and this is appealing mathematically. For the original Darboux transformation ([6]), the result is well 
known [19] 



Theorem 3 Denoting = (ipi,(f)i) T (i = 1,2,- •• N), then the N-fold Darboux transformation be- 
tween fields (|7|) can be reformulated as 

q[N] = q[0] - 2^, (16) 
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where 









-At (JV - 1} ^ •• 

1 ' 1 


■ -4 N ' 1] p N 


Ai = 


Y>1 


Aju ^0/v 


-ft 

a* (jv -% ■■ 

L T 1 


-0AT 




01 


0AT 


Pi 


Pn 




Af0i • 
Af- 2 ^i • 


A^0AT 

• A^"Viv 


xt N p 


*n N Pn 

x *(N-2) ,* 
./V ^Ar 


A 2 = 


ipi 


• A^~Vtv 


-0! 

A* (JV -% •• 






0i 


0AT 


Pi 


p n 



To find the determinant representations of our generalized Darboux transformation, we observe 
that the limit process presented in above subsection may be taken directly, namely (|10[) may be 
obtained by the following consideration 

q[2] = q[l} + lim 2(Q - C 2 )(P[2]) 2 i. 



Thus, as for the KdV case worked out in |19| . we may perform the limit on the determinant form (|16D 
and get 

Theorem 4 Assuming that N distinct solutions ^ = {ipi,cj)i) T (i = 1,2,- •• , n) are given for the 
spectral problem ([5a]) ([5b]) at £ = £i, • • • , C = Cn expanding 

(d + + 6) = CfVi + 1]<5 + 2]^ 2 + • • • + 0,[j, mi]8 mi +•••, 

(Ci + S) j (f>i(Ci + 6) = d<f>i + 0i[j", 1]<5 + 0,[j, 2]5 2 + • • • + 4 [j, mjr + • • • , 



with 



1 

m! 9£ r 



•[C%(C)]k c <: <k\iM 



i <9™ 

m! 9£ r 



■[C J 0i(C)]| 



fj = 0, 1, • • • , n, m = 1, 2, 3, • • • ), then we have 



D 



[N]=q-2-l, D 2 = det([H 1 ---H n ]), D 1 = det([Gi • • • G n ]), 



(17) 
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where N = n + Ylk=i m k an d 



Gi 



Hi 



c; 



N-l 



'h[N - l,mi 



4>i[0,mi 

(f>i[N, m, 
il>i[N-2,rm] 



Q N ^t 



*(iV-2) 



ipi[Q,mi] 
cj>i[N - l,m; 

0i [0, mj] 



Si 



C (iv-D 



-^i[JV-l,mi]* 

-^[O.m;]* 
^[JV-l.mi]* 

ijji[0,mi\* 
ijji[N,mi]* 
-</>i[N-2,mi]* 

-(pi [Q,mi\* 
iPi[N - l, mi )* 

ipi[0,mi]* 



We point out that, applying to special seed solutions, (fT7|) enables us to have a determinant form 
for higher order rouge wave solutions. We consider 



fa = i(C ie A - C 2 e~ A ), 4>i = (C 2 e A - C ie - A ), 



(18) 



where 



A = f^2Tp[x + \(i + f)t + Hf)], $(/) = EiIo^/ 2i > ^c. 



The associated Taylor expansions are 



p(i+fyM) = iVi(o) + ^ibM]r + --- + ^i[j,A r ]r v + ••• 



+ Wi(/)]|/=o, 



(2ra)!d/ 2n ' 

1 <9 2n 



(2n)! a/ 2 ™ 



[P(l + / 2 )^i(/)]|/=o, 



(j = 0, 1, • • • , iV, n = 1, 2, 3, • • • ). It follows that the iV-th order rogue wave solution for NLS equation 
([T]) reads as 



q[N] 



D 2 

1 - 2— 



(19) 
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where 



Di 



V>1 



i^-Vi 



^i[JV- 1,JV- 1] -(-i) 



(7V-1)^ 



7V-2 



■i/'i 



V>i[o,iv-i] 

<j>i[N -1,N -1] 
0l[O,iV-l] 

ifa[N-2,N- 1] 



^1 



-(-i) 



'I 



^i[0,JV-l] -<j>l 



l[0,JV-1] 



-0![7V-l,iV-l]* 

-^i[0,JV-l]* 
^i[JV-l,JV- 1]* 

^[0,JV-1]* 

^i[iV,iV-l]* 
-0i[iV-2,JV- 1]* 

-M0,N-1}* 
^i[JV-l,JV- 1]* 

^[0,JV-1]* 



For the case N = 2, above formula may provide the 2nd order rogue wave solution with two 
free parameters for the NLS equation, which has been analyzed in detail in [3]. It is shown that 
this solution splits into three 1st order rogue waves rather than two. Indeed, our results supply the 
high order rogue wave solutions with more free parameters, which determine the spatial-temporal 
structures of the solutions. In particular, the 3rd rogue wave solution, possessing four free parameters, 
may be worked out by setting N = 3 and $(/) = (b + ic)/ 2 + (e + ig)f 4 in our formula, whose explicit 
expression is omitted since it is rather cumbersome. In the following, we consider two special cases, 
which have different spatial-temporal patterns. 

• Case A 

In this case, we assume b = c = 0, |e| 3> or \g\ ^> 0. The corresponding 3rd order rogue wave 
solution is composed of six 1st order rogue waves, which array a regular pentagon. Interestingly, 
among the six 1st order rogue waves, one sits on the center and the rest locate on the vertices 
of the pentagon. After some calculations, we find that the radial distance from the center of the 
pentagon approximately equals to ^360 1//5 (e 2 + gr 2 ) 1 / 10 . For the case g = and e > 0, one of 
the vertices locates on the negative direction of x axis and the corresponding quadrantal angle 
for the (e, g)-plane is assumed to be zero. For the general non-zero g, the pentagon will rotate 
^9 along anticlockwise, where is the associated quadrantal angle for the (e,g)-plane. 



Case B 
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(a) The third-order rogue wave solution |g| 2 (b) Density plot for the third-order rogue wave 

solution \q\ 2 

Figure 1: The parameters e = 1000 and g = 

For the second case, we take the parameters e = g = 0, \b\ S> or |c| ^ 0. The correspond- 
ing third order rogue wave consists of the six first order rogue waves as well, which array an 
equilateral triangle. The distance between the center and any vertex approximately equals to 
90 1 / 6 (6 2 + c 2 ) 1 / 6 . For the case c = and 6^0, one of the vertices locates on the positive 
direction of x axis and and the corresponding quadrantal angle for the (6, c)-plane is assumed 
to be zero. For the general non-zero c, the triangle will rotate ^9 anticlockwisely, where is the 
related quadrantal angle for the (b, c)-plane. When b = c = e = g = 0, the rogue wave is the one 
considered by Akhmediev et al in [2]. 

Let (xi,ti) (i = 1, 6) be the coordinates of the six peaks, then we find that the third order rogue 
wave may be approximated by 

e l \ di = 1 + 4(> - Xi f + 4(t - U) 2 , 

the "center of mass" is at the origin (0, 0) in both cases. With raising of the order, the rogue solution 
contains more free parameters and exhibits more interesting spatial-temporal structures. For instance, 
choosing the proper parameters, the 4th order rogue wave possesses the regular heptagon spatial- 
temporal pattern (Fig. [3]), and a second order rogue wave locates on the center. Naturally, we 
conjecture that the spatial-temporal pattern of N-th order rogue wave possesses 2N — 1-gon spatial- 



q[3] * £ 



-1 + 4 



1 + 2i(t - U 

di 
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(a) The third-order rogue wave solution |g| 2 (b) Density plot for the third-order rogue wave 

solution \q\ 2 

Figure 2: The parameters b = 100 and c = 

temporal pattern by choosing the parameters $*(/) = cf 2i - N ~ 1 ^ (|c| 3> 0) in formula (I19p . where c is a 
complex number. 

We conclude this section with the following remarks in order: 

1. If we choose the seed solution ^[0] = 0, we may obtain the higher soliton solution |1 1| 121] , 

2. The integrable Hirota equation 

i<?t + T^Qxx + \q\ 2 q - ia(q xxx + 6\q\ 2 q x ) = 0, a is real constant (20) 

is the third flow of NLS hierarchy. Its rogue wave solutions and rational solutions were discussed 
in pQ. To obtain the iV-th order rogue wave solution for Hirota equation (j20|) from ()19j) . we 
merely need to modify the A = fyj2 + f 2 [x + i(l + f 2 )t + a{2 + 4(1 + f 2 ) 2 )t + $(/)] in (fl8j) . 
Rather than giving the explicit exprssions, which is lengthy we plot the third-order rogue wave 
solution for Hirota equation with the Fig. HI Due to the third-order dispersion and time-delay 
correction to the cubic term [1], the solution enjoys the different behavior from NLS equation. 

4 Conclusion 

Through a limit procedure, we have generalized the original Darboux transformation for the NLS 
equation. This Darboux transformation, in particular, allows us to calculate higher order rogue wave 
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(a) The fourth-order rogue wave solution \q\ 2 (b) Density plot for the fourth-order rogue wave 

solution | q \ 2 

Figure 3: The parameters $(/) = 5 x 10 5 i/ 6 




solutions in a unified way. We believe that the idea is rather general and could be applied to other 
physically interested models as well. These spatial-temporal structure for N-th order rogue wave may 
be useful to research spatial-temporal distributing of rogue wave in deep water. 
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